The Hypertriton in Effective Field Theory 



H.-W. Hammeifl 

Department of Physics, The Ohio State University, Columbus, OH 43210, USA 

(Dated: October 15, 2001) 

Abstract 

Doublet Ad scattering and the hypertriton are studied in the framework of an effective field 
theory for large scattering lengths. As in the triton case, consistent renormalization requires a one- 
parameter three-body force at leading order whose renormalization group evolution is governed by 
a limit cycle. Constraining unknown parameters from symmetry considerations and the measured 
binding energy of the hypertriton, we calculate the low-energy phase shifts for doublet Ad scatter- 
ing. For the low-energy parameters, we find = (16.8124) ™d rAd = (2.3 ± 0.3) fm, where 
the errors are due to the uncertainty in the hypertriton binding energy. Since the hypertriton is 
extremely shallow, low-energy three-body observables in this channel are very insensitive to the 
exact values of the AA^ low-energy parameters. 
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I. INTRODUCTION 



In recent years, there has been much interest in applying Effective Field Theory (EFT) 
methods to nuclear systems with two or more nucleons |jl|, H ^, EFT's provide a 
powerful framework to explore a separation of scales in physical systems in order to perform 
systematic, model- independent calculations 0. If, for example, the momenta k of two 
particles are much smaller than the inverse range of their interaction 1/R, observables can 
be expanded in powers of kR. All short- distance effects are systematically absorbed into 
a few low-energy constants using renormalization. The EFT approach allows for accurate 
calculation of low-energy processes with well-defined error estimates. For low-energy nuclear 
few-body systems, the long-distance scale is set by the large two-body scattering lengths, 
while the short distance scale is set by the range of the nuclear force or the inverse pion 
mass. In an EFT, the long-distance physics is included explicitly, while the corrections from 
short- distance physics are calculated in an expansion in the ratio of these two scales. 

For very low energies and momenta (p<m^), even pion exchange can be considered 
"short-distance" physics. In this case, one can use an effective Lagrangian including only 
contact interactions. The large S-wave scattering lengths require that the leading two-body 
contact interaction is treated nonperturbatively [0, § . In the two-body system, this program 
has been very successful (see e.g. Refs. ||, |10| and references therein). 

In the nuclear three-body system, considerable progress has been made as well [pT| , pT^ , 
T3|, O, O, O, O, EOi . Most of the work has concentrated on the neutron-deuteron 



system. In channels where the Pauli principle or centrifugal barrier suppress sensitivity to 
short- distance physics, the EFT can be extended in a straightforward way |Tl|, O, Il3|, [14 



Recently, Coulomb repulsion has been included and the phase shifts for proton-deuteron 
scattering in the spin quartet (J = 3/2) have been calculated [0. Most interesting, however, 
is the S-wave in doublet {J = 1/2) neutron-deuteron scattering which displays a number of 
surprising phenomena such as the Thomas and Efimov effects or the Phillips line ^ ^ 



In this channel, the renormalization requires a one-parameter three-body force at leading 
order whose renormalization group evolution is governed by a limit cycle [|T6], ^ . Variation 



of this three-body force gives a compelling explanation of the Phillips line P2[. Recently, 



the effective range corrections to S-wave neutron-deuteron scattering in the doublet channel 
have been calculated At this order, good agreement with available phase shift analyses 
and potential model calculations is obtained. 

The purpose of the present paper is to study doublet Ad scattering and the hypertriton 
using EFT methods. The hypertriton is the lightest hypernucleus and consists of a triton 
with one neutron replaced by a A. It is the simplest strange three-body halo nucleus with a 
separation energy into a deuteron and a A of only = 0.13 ± 0.05 MeV [0, |25|. The total 
binding energy is = 2.35 MeV. The hypertriton has been studied extensively using vari- 
ous hyperon-nucleon potentials (see e.g. Refs. [|^, ^7\, |2^, |3l| and references therein). 
The most sophisticated calculations include both tensor forces and A ^ S conversion effects 
27| . Obtaining the correct binding energy for the hypertriton in these potential models is 



delicate, which has motivated the use of this system to study details of the hyperon-nucleon 
interaction |]30 |. 

The EFT framework offers a different perspective and stresses the universal aspects of the 
problem. In this paper, we study Ad scattering and the hypertriton in the framework of an 
EFT for three-body systems with large two-body scattering lengths [|1^, |1^]. The small A^— A 
mass difference can be exploited by expanding in the parameter y = (Ma — M) / (Ma + M) ^ 
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0.1. After this expansion is performed, the hypertriton emerges naturally as a shallow Efimov 
state. To leading order, both NN and AA^ S-waves will contribute to the hypertriton. While 
the A^A^ S-waves have a large scattering length, the AA^ scattering lengths are believed to 
be natural. Therefore, it is not obvious that an EFT for large two-body scattering lengths 
is applicable. In the following, we motivate the use of this EFT. 

Unfortunately, there are only few AA^ scattering data and these are at relatively high 
energies. Model independent analyses using effective range theory indicate mainly S-wave 
scattering, but afford many solutions for the low-energy parameters and are essentially 
inconclusive . For example in Ref . |Q , the extracted one-standard-deviation bounds 



were > ai > — 15 fm, < ri < 15 fm, —0.6 fm > as > —3.2 fm, and 2.5 fm < ra < 15 
fm where ai (as) and ri (rs) are the singlet (triplet) scattering length and effective range, 
respectively. These bounds do not exclude the case of large scattering lengths. 

Various sets of low-energy parameters are also known from extrapolations of the few 
higher energy data using hyperon-nucleon potentials ^ |36|, ^ . The extracted AA^ 



scattering lengths and effective ranges are generally of natural size [0(1/171^^)]. Both the 
^5*1 and ^5*0 AA^ partial waves have a virtual bound state. The pole momenta extracted 
from the low-energy extrapolations are ~ 60 MeV for the ^5*0 and ~ 70 MeV for the ^5*1 
partial wave. However, it is not clear how large the errors and model-dependence in these 
low-energy extrapolations are. 

An EFT analysis of hyperon nucleon scattering and hyperon mass shifts in the nuclear 
medium was recently performed in Ref. [^]. The AA^ scattering lengths, however, were used 



as input in this analysis. 

Because the hypertriton is so weakly bound, it should only be sensitive to the extreme 
long-distance properties of the AA^ interaction. The typical momentum scale of the hyper- 
triton can be estimated from its binding energy via 



7s^ ~ 2^{MB^ - jj)/3 ^ 0.37. , (1) 

where M = 938.9 MeV is the nucleon mass and 7. = 45.68 MeV the deuteron pole momen- 
tum. Since 7;^ is so small, the contribution from the effective ranges will be suppressed even 
if the AA^ low-energy parameters are 0{l/mT^) as extracted from the potential models (see 
Section |^ for a more quantitative estimate). Consequently, the use of the EFT for large 
two-body scattering lengths is justified. 

In the next Section, we write down the effective Lagrangian and discuss the power count- 
ing of the EFT. We then review the EFT for the A^A^ and AA^ two-body subsystems. In 
Section |T|, we obtain the three-body integral equations for the ANN system and in Sec- 
tion 1^ we present our results and conclusions. Some technical details are given in the 
Appendices. 



II. TWO-BODY SYSTEM 

The power counting in an EFT is determined by the physical scales in the system under 
consideration. For the hypertriton, we are only interested in very low-energies, where all 
physics (even pion exchange) can be considered "short-range". The two-body subsystems 
(A^A^ and AA^) are characterized by two scales: a long-distance scale (the large S-wave scat- 
tering length) and a short-distance scale that is given by the longest-range physics excluded 
from the EFT. For an EFT without pions, one expects the short-distance scale R to be of 
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C(l/m7r). However, since the A (J = 0) and nucleon (J = 1/2) cannot exchange a pion 
with 1=1 and conserve isospin, the short-distance scale for the hypertriton should be set 
by the two-pion exchange |]2^. A power counting to deal with unnaturally large two-body 



scattering lengths has been proposed in Refs. 0, §]. This counting takes p ~ 1/a ~ Q, where 
p is the typical momentum and 1/a is the inverse scattering length. The pole or binding 
momentum 7 is of 0{Q) as well since 7 = 1/a -|- (9(r/a^) with r R the effective range. 
The expansion of the EFT is in powers of QR ~ 7r (see Refs. ^ for more details). 
For the study of three-body systems, it is convenient to employ the dibaryon formalism 
where an auxiliary field is used to describe two baryons interacting in a given partial 



wave. We write down an effective Lagrangian including nucleons, A's, and dibaryon fields 
for the deuteron {d) as well as the ^5*1 {v?) and ^Sq {u^) AN partial waves, 

/:^iVt(.c^. + ^)iV + At(.5, + ^jA (2) 
+Aa4di - ^ [djN^iiT2)itaia2)N + H.c." 

+Ai{uyu^ - igi [{uyA^{ia2)N + H.c] + . . . , 

where H.c. denotes the Hermitian conjugate and the dots indicate terms with more deriva- 
tives and/or fields. The terms with more derivatives are suppressed at low energies, while 
four- and higher-body forces do not contribute. Three-body terms will be addressed in 
the next Section. In Eq. (^, indices occuring repeatedly are summed over and cTj (xj) are 
the usual Pauli matrices acting in spin (isospin) space. Spinor and isospinor indices are 
suppressed. The parameters A and g are not independent and only the combination A/g"^ 
enters physical observables. 

The Lagrangian (^ describes low-energy S-wave scattering of nucleons and A's and re- 
produces the effective range expansion up to the scattering length term, which is sufficient 
for our purposes. To go to higher orders in QR, the two-body effective range can be included 
with a kinetic term for the auxiliary field [^. Only the partial waves contributing in the 



hypertriton channel to leading order in QR are included in the effective Lagrangian (^. 
The ^5*0 A^A^ amplitude does not contribute because of isospin, while higher partial waves 
and the tensor force are suppressed by at least two powers of the expansion parameter QR 
1^. After integrating out the auxiliary fields d, u^, and m^, it is straightforward to show 



that the Lagrangian (^ is equivalent to a Lagrangian with nucleon and A fields only (cf. 
Refs. 111,113). 

Since the theory is nonrelativistic, all particles propagate forward in time and the nucleon 
and A tadpoles vanish. The propagator for nucleons and A's is simply 

iS{po,p) = . , (3) 

Po — P i (,2mj + le 

where m is the mass of the corresponding particle. The dibaryon propagators are more 
complicated because of the coupling to two-baryon states. The bare dibaryon propagators 
are constant, i/A, but the full propagators are dressed by baryon loops to all orders, as 
illustrated for the AA^ state in Fig. |l|. All diagrams in Fig. |l| are of the same order because 
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FIG. 1: Dressing of the bare AA^ propagator. The dashed hne indicates a A, the sohd hne a 
nucleon, the thick sohd hne a bare, and the double line a full dibaryon propagator. 

2'n'A/{Mg'^) = 1/a ~ Q. Summing the resulting geometric series leads to 

/ ^ 27r —i 

iDj{po,p) = 2 / , (4) 

t'^a] -7, + {po - pV [2(M + Ma)]) - te 

where /xa = MM/^/ (M + Ma) is the reduced mass of the AN system, j G {1, 3}, and 7^ is 
the corresponding pole momentum. Divergent loop integrals are regulated using dimensional 
regularization. The full propagator for the deuteron is [|12| , |1^ 



/ ^ 27r —i 
iDd{po,p} = / , (5) 



where 7^ = 45.68 MeV is the deuteron pole momentum. The values of 71 and 73 will 
be discussed later, together with the three-body results. The scattering amplitudes are 
obtained by attaching external baryon lines to the full propagators from Eqs. (§J^) [|13] • AH 
dependence on the bare coupling constants ga, gi, and g^ cancels in observable quantities. 



III. THREE-BODY SYSTEM 



We now apply the Lagrangian (0) to the ANN system in the J = 1/2 channel, which 
has the hypertriton as a three-body bound state. We do not include explicit A S con- 
version effects Such effects are accounted for by a ANN three-body force to be in- 
troduced in the following. Explicit S degrees of freedom have been integrated out from 
the effective Lagrangian. This procedure is justified because the characteristic momentum 



7^ < VMa(Ms - Ma) ^ 300 MeV |0| 



To leading order in QR, only relative S-waves contribute. A ANN state with the quantum 
numbers of the hypertriton ( = ^~^) can then be constructed in three different ways: (i) a 
^5*1 AN partial wave plus another nucleon, (ii) a ^Sq AN partial wave plus another nucleon, 
and {in) a ^Si NN partial wave (a deuteron) plus a A. The ^Sq NN partial wave does 
not contribute because of isospin. The leading correction comes from the two-body effective 
ranges which enter at 0{QR). Contributions from higher partial waves and the tensor force 
are suppressed by at least two powers of QR p^ , pO| ; in potential models, their contribution 



is found to be small as well 

As a consequence, we need three three-body amplitudes, Ta, Tb, and Tc, to describe the 
hypertriton. The amplitudes satisfy the coupled integral equations shown in Fig. ^. The 
amplitude Ta describes J = 1/2 Ad scattering, while the remaining two amplitudes have a 



See Refs. [Eol ETl for a discussion of a similar issue in the A system. 
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FIG. 2: Coupled integral equations for J = 1/2 Ad scattering. Nucleons (A's) are indicated by 
solid (dashed) lines, respectively. The d dibaryon is indicated by double solid lines, while the 
(^S'o) AA^ dibaryon is indicated by the dashed-solid double line with index 3 (1). 



Ad initial state but go into the final states (i) and (ii) from above. Note also that there is 
no tree level diagram for Ta because strangeness is conserved. 

The details of the derivation of the integral equations are given in Appendix At a 
first glance, it appears that the Efimov physics known from the neutron-deuteron system 
||16|, llTl |23| is not present here. The masses of the A and nucleon are different and the 



integral equations ( [A8| ) are not scale invariant for large loop momenta. However, since the 
A — N mass difference [characterized by the ratio y = {M\ — M)/{Mj^ + M) ^ 0.1] is 
small, the equation becomes approximately scale invariant. Therefore, the Efimov physics 
will be present and can be made manifest by expanding around the limit y = 0. The error 
introduced by keeping only the leading order in this expansion is at most of the order of 
the effective range corrections, which we neglect. Corrections to this limit can be calculated 
within the EFT. Setting y = in Eq. (|A8D, we obtain: 



1 M 



TT JO 

M 

1 M 



dq {Lb{p, q, E)TB{k, q) - 3Lc(p, g, E)Tc{k, q)} 



(6) 



Tcik,p) 



TT Jo 



1 /"^c 

Lj{p, k,E) / q^ dqlAip, q, E)TA{k, q) 

TV JO 

q^ dq {Lb{p, q, E)TB{k, q) + 3Lc(p, g, E)Tc{k, q)} 



1 r^'= 

Li{p, k,E) + - q^ dqlAip, g, E)TA{k, q) 
M vr Jo 



TT JO 



Ac 



q^ dq {Lb{p, q, E)TB{k, q) - Ldp, q, E)Tc{k, q)} 



where k {p) denote the incoming (outgoing) momenta in the center-of-mass frame and A^ is 
a momentum cutoff introduced to regulate the integral equations. In Eqs. (|), the outgoing 
momentum p is taken off-shell. The total energy is E = (4M) — 7J/M. 
In the limit y = 0, the functions Lj and La are given by 



1 f e+pk + p'-ME \ 



(7) 
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1 (q^+pq+£-ME\ 

La{p, q,E) = — In — — — - 

pq \q'^ — pq + — Mtj J 




ME 



It 



Lb (Lc) are identical to La with 7^ replaced by 73 (71), respectively. The amplitude Ta{p, k) 
is normalized such that 

Stt 1 

with 5 the elastic Kd scattering phase shift. 

The integral equations (^) show the same feature as in the case of J = 1/2 neutron- 



deuteron-scattering or three spinless bosons [0, |T^. Without the cutoff A^, their solution 
would not be unique The origin of this nonuniqueness can be understood by solving 
Eq. (|^) for asymptotically large off-shell momenta ^ g ^ 7^ ~ A;. In the asymptotic 
limit, the inhomogeneous terms can be neglected and the equations become independent of 



Ta{p) 
Tb{p) 
fcip) 




p2 _ pq q2 



{fBiq)-3fciq)} 
[fAiq) + fBiq) + 3fc{q)} 
\{fA{q)-fB{q)+fc{q)} , 



(9) 



where we have defined Tj{p) = pTj{k ~ 7d,p) for j G {A, B, C}. Note that the asymptotic 
equations do not determine the overall normalization of the amplitudes, which will be given 
from matching to the low-energy solution of the full equation. The following discussion, 
however, is independent of this normalization. We can decouple the equations by defining a 
new set of amplitudes Ti, T2, and T3 via 



(10) 
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-I-V5 -1 + V5\ 
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The new amplitudes Ti, T2, and T3 fulfill the equations 
Ti{p) = 
T2{p) = 

Up) = 




(11) 



Tsiq). 



All three equations are scale invariant and symmetric under the inversion q 
Ref. [H2|, it was shown that an equation of the type 



fip) 



4A 
\/37r Jo 



In 



p + pq + q 

p2 _ pg _|_ g2 



1/q. In 
(12) 
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has a unique solution in the form of a power law if A < Ac = 3-\/3/(47r) ^ 0.4135 (see 
also Refs. |]TB|, [T^). This condition is clearly satisfied for the first two equations. The 
equation for T3, however, has A = (1 + a/5)/2 ^ 1.62. As a consequence, there are two 
linearly independent complex solutions T^lkyp) = where si = 1.35322. The relative 

phase between these two solutions is not determined by Eq. (0). For a finite cutoff Ac this 
phase is fixed but strongly depends on Ac. This Ac dependence can absorbed by adding a 
one-parameter three-body force in the equation for T3, 

TM = ^i±^ ^ {,„ f Pl±Pl±l) + 2HiAM , (13) 

x/37r 2 Jo Q \ - PQ + A. 



where H(Ac) is a dimensionless function of the cutoff. In Refs. [T^, |T3|, it was shown that 
if the three-body force runs with cutoff as 

fr(.\) - ^^^^^^ H^c/K) - arctg(l/ gQ) 
^ sin(siln(Ac/A,) + arctg(l/si)) ' ^ ^ 

all low-energy observables are independent of Ac. The renormalization group evolution of 
H{Ac) is governed by a limit cycle. In particular, if the cutoff is increased by a factor of 
exp(7r/si) ~ 10.2, H{Ac) returns to its original value. A* is a dimensionful parameter that 
determines the asymptotic phase of the off-shell amplitude [|16|, |1^. It has to be fixed from 



a three-body observable. Once its value is known, all other observables can be predicted.^ 
Note, however, that A^, is only determined up to factors of exp(7r/si). 

Formally, the three-body force term is obtained by adding a nonderivative three-body 
contact term to the effective Lagrangian (0). The form of this term can be determined 
by starting from a general structure and matching the unknown coefficients to reproduce 
Eq. ([TB|). The exact expression for this three-body term is given in Appendix It is 
interesting to note that the value si ~ 1.35 differs from sq ~ 1.006 which characterizes the 



three-nucleon force in the triton channel ^7j. This does not lead to a contradiction because 
the hypertriton channel (/ = 0) and the triton channel (1 = 1/2) have different isospin. 

For the numerical studies, it is useful to write a renormalized integral equation by choosing 
a special cutoff. 



An = A^exp 



1 / /I 
— nvr + arctan — 

Si V Vsi 



(15) 



with n > an integer, at which the the three-body term proportional to H in Eq. (|T3D 



vanishes [|1^]. This procedure does not remove the dependence on the three-body parameter 
A* but rather moves its dependence into the cutoff. All observables are independent of n, 
but it should be chosen such that A„ ^ and finite cutoff corrections of order jdZ-^n can 
be neglected. In the following, we will be using Eqs. (^ with Ac = Ai. The Ad scattering 
amplitude is obtained by numerically solving Eqs. @ for the desired value of the total 
energy. The hypertriton binding energy can be extracted by solving Eqs. (|^) for negative 
energies and locating the position of the bound state pole in T4. 



^ One can either specify the dimensionless couphng H at a specific cutoff Ac or the dimensionful low-energy 
parameter A*. This is similar to dimensional transmutation in QCD. 
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FIG. 3: Values of A-^^/^d and 7j/7d that lead to the correct binding energy for the hypertriton. 



IV. RESULTS 



In order to apply the EFT to the ANN system, we have to fix the parameters in Eqs. (^. 
The deuteron binding momentum 7^ = 45.68 MeV is known very well. As mentioned above, 
the low-energy parameters for the AA^ S-waves cannot be extracted unambiguously from a 
model-independent analysis of low-energy scattering data P^ , It is known, however, 
that the AA^ system is unbound, which requires 71 and 73 to be negative. The errors in 
the low-energy parameters extracted from extrapolations using hyperon-nucleon interaction 
potentials ^ are not well known and can only be estimated from the variation 

with different potentials. Therefore, we do not take the extracted parameters as input. 
Rather, we vary the AA^ pole momenta over the range of validity of the EFT and study the 
implications on the hypertriton and Ad scattering. 

In principle, we have three unknown parameters: 71, 73, and the three-body parameter 
A*. However, we will first address the simpler scenario 71 = 73, which follows from the SU{6) 
spin-flavor symmetry of the quark model |^3|, Such an SU{6) symmetry also emerges in 



the large- Ac limit of QCD Symmetry violating corrections are suppressed by 1/Nc and 
the strange quark mass. Since both isospin and strangeness of the ^5*1 and ^5*0 AA^ partial 
waves are the same, the difference between 71 and 73 is a pure spin splitting effect. In most 
hyperon-nucleon potentials this effect is of the order 30% ||3 
splitting leaves us with two parameters: 7j(= 71 = 73) and A^ 



55], |3|, 0. Neglecting this 
In the following, we vary 
these two parameters and study the constraints that follow from a correct description of the 
hypertriton. We will also elucidate the consequences for Ad scattering. 

First, we study the constraints on the parameters 7^ and A* by the requirement that 
the hypertriton binding energy be reproduced. An interesting question is whether there are 
certain values of 7^ and A,, that are incompatible with the hypertriton data. Furthermore, 
it is possible that the approximation 71 = 73 proves too restrictive. In Fig. |^, we plot A^, 
as a function of 7^ with the requirement that the hypertriton binding energy is reproduced. 
We vary 7j from to —2'-/^, which is about where corrections of 0{QR) are expected to 
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FIG. 4: Ad scattering phase shifts for different values of 7^ and that lead to the correct hyper- 
triton binding energy (all quantities in units of 7^). 



become important. Positive values of 7^ are excluded because the AA^ system is not bound. 
All allowed values fall on a curve with a linear dependence on 7^, 



A,/7, = 1.08 - 0.557i/7d. 



(16) 



For every value of 7^ in the interval a corresponding value of A* that reproduces the hyper- 
triton can be found. Clearly, the exact value of 7^ and therefore the AA^ scattering lengths 
cannot be determined from the hypertriton binding energy. This is in contrast to Ref. pO 



where it was found that with standard hyperon-nucleon potentials, the correct binding en- 
ergy for the hypertriton cannot be reproduced if the singlet AA^ scattering length differs 
more than 10% from 1.85 fm. 

In Fig. ^ we compare the scattering phase shifts for Ad scattering in the hypertriton 
channel for different points on the curve in Fig. ^. Fig. ^ shows three curves with 7j chosen 
to be equal in magnitude to 7^ as well as a factor of two smaller and larger. The phase shifts 
agree very well for small momenta around threshold and begin to differ slightly close to 
the deuteron breakup threshold. This implies that low-energy Ad scattering does not give 
additional constraints for the values of 7^ and A*. Consequently, neither the hypertriton 
binding energy nor the low-energy Ad phase shifts can be used to fix the exact value of 7^. 
Such information can only come from low-energy AA^ scattering. On the other hand, this 
observation allows for a prediction of the low-energy Ad phase shifts [given by the curves in 
Fig. 1^ independent of the exact values for 7j and A,,. This is a direct consequence of the 
shallowness of the hypertriton. Extracting the low-energy doublet Ad scattering parameters 
from the phase shifts, we find for the scattering length a^d and effective range r^d'- 



= (16.81^:^) fm and r^a = (2.3 ± 0.3) fm , 



(17) 



where the errors are due to the uncertainty in the hypertriton binding energy. Our value for 
the scattering length is in good agreement with Ref. while the effective range is about 
1 fm smaller. 
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FIG. 5: Phillips line for different values of 7^ (all quantities in units of 7^). 



It is also interesting to look at the "Phillips line" , which describes the correlation between 
the Ad scattering length a\d. and the hypertriton binding energy B^. The Phillips line was 
observed in the neutron-deuteron system, where the results of calculations of the doublet 
nd scattering length and the triton binding energy using different A^A^ potentials fall on a 
line when plotted in a plane A similar line exists for the hypertriton (cf. Fig. 7 in 

Ref. [^). In the EFT framework, the Phillips line is parametrized by variations in the 
three-body parameter A* [0, O. In Fig. ^ we show the Phillips line for the three values 
of 7j from Fig. ^. For small B^, the different Phillips lines coincide exactly (the physical 
hypertriton corresponds to MB^ ^ 1.067^) and deviate from each other only at extremely 
large binding energies. However, the EFT is expected to break down when the three-body 
binding momentum 7^ from Eq. (0) becomes of the order of the pion mass. This is the case 
for MB^ > 87J. Consequently, the region where the Phillips lines differ is unphysical. For 
all practical purposes the Phillips line is independent of 7^. 

The characteristic momentum that sets the scale of corrections to the large-scattering- 
length approximation is given by the binding momentum of the hypertriton. Using the 
measured value for B^ and Eq. (|l|), we obtain 7^^ ^ 0.3 7^. This momentum can be used 
to estimate the contribution of the effective range term in the two-body amplitudes which 
is the leading correction at 0{QR). For our approximation to be accurate, we need I7I ^ 
|r(7^)^/2| in each of the two-body partial waves. Since (7^)^/2 pa 0.05 7J is extremely small, 
the effective range term gives only a small correction even when the scattering length and 
effective range are of the same order. 

We have also studied the more general case 71 7^ 73. The conclusions in this case are 
unchanged: requiring the correct binding energy of the hypertriton does not constrain the 
specific values of 71 and 73 as long as they are within the range of validity of the EFT 
(0 > 7i, 73 > — 27rf). The results for Ad scattering confirm Fig. | and Eq. (|T^. 

The insensitivity to the precise values of 71 and 73 and the small separation energy of 
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the hypertriton into a deuteron and a A {B^ = 0.13 MeV) can be exploited even further: it 
is possible to write an EFT with only deuterons and A's as degrees of freedom |4^. Effects 
from the deuteron structure and AA^ interaction can be included perturbatively via local 
operators. 

To summarize, we have discussed the hypertriton and doublet Ad scattering in an EFT 
for large two-body scattering lengths. Due to the small binding energy of the hypertriton, 
our analysis is expected to hold even if the AA^ low-energy parameters are (9(l/m^) as 
extracted from potential models. As in the triton case, consistent renormalization requires 
a one-parameter three-body force at leading order whose renormalization group evolution 
is governed by a limit cycle. The period of the limit cycle exp(7r/si) ^ 10.2 is a factor 
two smaller than for the triton. In contrast to standard potential models where the binding 
energy of the hypertriton is very sensitive to the singlet AA^ scattering length ||3^, the three- 
body force can always compensate for changes in the low-energy AA^ parameters. While the 
AA^ low-energy parameters cannot be fixed by requiring the correct binding energy for the 
hypertriton, the shallowness of the hypertriton still allows for a unique prediction of low- 
energy doublet Ad scattering. Using the experimental binding energy of the hypertriton 
as input, we find for the scattering length a^d = (16.8l2;4) and for the effective range 
= (2.3 ± 0.3) fm. In the future, it would be interesting to go to higher orders and apply 
EFT methods to a wider class of observables and other halo nuclei. Work in this direction 



IS m progress 
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APPENDIX A: DERIVATION OF THREE-BODY EQUATION 

In this Appendix, we give some details of the derivation of the integral equation (P) for 
the general case Ma 7^ M. From the Lagrangian and the Feynman diagrams in Fig. |^, 
we obtain^ 

m p)., = g,g, J ^ t^{k, g)., /fJ^'J'f'/^'^. (Al) 

^ ^ ^ 2Ma 2M 2M ^~ 

^ ^ 2Ma 2M 2M ~^ ^ 



^ 2Ma 2M 2M ''^ 



V / ^ O J\yf OA/f. Oil//' ~r ''^ 



2M 2Ma 2M 



See Refs. ||T^, for details on the evaluation of the Feynman diagrams. 
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^ ' ^ 2M 2Ma ^ 

^ ^ 2M 2Ma ^ 



2Ma 2M 2M 



^ 2M 2Ma 2M ~ ''^ 



^ 2M 2Ma ^ 

^ ' ^ 2M 2Ma 

where i,j,k are vector indices and Ci,/3,^ are spinor indices in spin space. k{p) are the 
incoming (outgoing) momenta in the center-of-mass frame and E — k'^/[AM) — 7J/M + 
/c^/(2Ma) is the total energy. The momentum p is taken off-sheU. The isospin dependence 
of the amphtudes Ib and tc has been absorbed into the amphtudes via the definition 

t%{Kp).p = 4ik,P}U^2)be, (A4) 

and similarly for tc, where e, b are isospinor indices. Next we project onto total angular 
momentum J = 1/2 by using the definitions 

tA/Bik,p) Sa/3 = ^{cri)aa'fi^/B0,P)a'P'{crj)/3rp (A5) 

and projecting on relative S- waves. We also have to account for the wave function renor- 
malization of the deuteron, 

(iDJpo,p))~^ , = -. (A6) 



Defining 



TA{k,p) = ZatA{k,p) (A7) 
TBik,p) = Za^tB{k,p) 

93 

Tc{k,p) = Zd—tc{k,p) , 
91 

we finally obtain the integral equations 

Ta{Kp) = - ^^i\y^ l"^^ dq {Lb{p, g, E)Tb{K q) - ^Lc{p, g, E)Tc{K q)} (A8) 
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Tc{k,p) 



M 7T Jo 

1 r^'' 

71 7 / 9^ dq {Lsip, q, E)TB{k, q) + 3Lc(p, g, E)Tc{k, q)} 

7r(l - y) Jo 

^^Li{p, k,E) + - f q^ dqLAip, q, E)TA{k, q) 
M TX Jo 



1 



7r(l - y) Jo 



Ac 



rfg {Lb{p, q, E)TB{k, q) - Lc{p, q, E)Tc{k, q)} 



where Ac is a momentum cutoff and y = (Ma — M) / (Ma + M) ^0.1 characterizes the A — 
mass difference. The functions Lj, La, Lb, and Lb are given by 



Lrip,k,E) 
LA{p,q,E) 



1 ^^fk'^/{l + y)+pk+p'^-ME\ 
^ ^ \F/(1 + y)-pk+p^ - ME J 
q^ / {I + y) + pq + p^ - ME \ 
pq \q^ /{I + y) — pq + p"^ — ME j 



X 



-Id + (3 - y)/{l + y)/4 - ME - 



LB{p,q,E) 



]_^^( q^+p^+ pq{l -y)- ME{1 + |/) ' 
\q^+p^— pq{l ~ y) — ME{1 + y) 



X 



-73 + \/g2 (3 + 2|/ - 1/2) /4 - ME(1 + y)-ie 



■1 -1 



LB{p,q,E) 



]_^^( q^ + pq + p^ / {I + y) - ME \ 
pq \q'^—pq + p^/{l + y) — MEj 



X 



73 + V?' (3 + 22/ - 2/2) /4 - ME{1 + y)-ie 



Lc [Lc) is identical to Lb (Lb) with 73 replaced by 71. 



APPENDIX B: THREE-BODY FORCE TERM 



(A9) 



The three-body term proportional to H(Ac) in Eq. (|T^ can be obtained by adding a 
local, nonderivative three-body contact term to the Lagrangian (^). This term is most 
easily determined by starting from a general structure and matching the coefficients of the 
individual terms to reproduce the three-body term in Eq. (|13D . This leads to, 



MH{A^ 



lQV5gl kA 



+ 2(5 + V5)gdg3 [d*iNaai(ri(ri')af3A*f3ir2)abiul)b + H.c. 
- (5 + 3V5)gig3 [iuf)lN^a{ai)af3N;,iu')b + H.c." 

+ 2(5 + V5)gdgi [diNo,a{(Ti)apAl{T2)abiu^)b + H.c. 



(Bl) 
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where a, P (a, b) are spinor (isospinor) indices and /, /' are vector indices. Eq. (|B1|) represents 
a contact three-body force written in terms of dibaryon, nucleon, and A fields but to leading 
order is equivalent to a three-baryon contact force. This becomes obvious if the dibaryon 
fields are integrated out by performing a Gaussian path integral ill, Il6|, IT? . 
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